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BIOCHEMICAL AND MORPHOMETRIC 
RELATIONSHIPS AMONG SOME 
MEMBERS OF THE CARDINALINAE 

Jeffrey W. Tamplin, 2 James W. Demastes, 1 - 2 and 

J. V. Remsen, Jr 1 2 

Abstract. —We used starch-gel electrophoresis to analyze relationships among 1 3 species 
in the subfamily Cardinalinae. These results were compared to a morphometric analysis of 
15 skeletal characters and to a previous morphometric analysis of the subfamily by Hellack 
and Schnell (1977). Our results supported the phenetic classification of Hellack and Schnell 
and were more consistent with the classification of Hellmayr (1938) than with that of Paynter 
(1970). Divergence in morphometries of skeletal characters in the Cardinalinae has been 
primarily in size. Electrophoretic data suggest that the Yellow-shouldered Grosbeak ( Cary- 
othraustes humeralis) is not a member of the Cardinalinae. In addition, if the broad genus 
Saltator is to be retained, then the genus Pitylus should be merged into it, and the Dickcissel 
(Spiza americana ) is the outgroup to all other Cardinalinae examined (excluding Cary- 
othraustes humeralis). Received 22 Nov. 1991, accepted 15 April 1992. 


The subfamily Cardinalinae (Emberizidae/Fringillidae), consists of 37- 
42 species of cardinals, grosbeaks, and buntings (Paynter 1970, Morony 
et al. 1975, Sibley and Monroe 1990). Following a trend over the last 
three decades toward broader generic limits in avian systematics, seven 
genera (Hedymeles, Richmondena, Pyrrhuloxia, Cyanocompsa, Cyano- 
loxia, Guiraca, and Porphyrospiza ) recognized by Hellmayr (1938) were 
merged into other genera by Paynter (1970). However, as has been typical 
of avian systematics throughout much of its history, reasons for taxonomic 
changes were not explicit. 

Hellack and Schnell (1977) used plumage characters and skeletal mea¬ 
surements from Hellack (1976) to investigate phenetic relationships among 
the Cardinalinae. Although their phenetic classification resembles those 
of Hellmayr (1938), Paynter (1970), and Sibley and Monroe (1990), Hel¬ 
lack and Schnell’s groupings were generally more similar to those of Hell¬ 
mayr and of Sibley and Monroe than to those of Paynter. For example, 
morphometric results suggested that Paynter’s expanded genus Passcrina 
was a heterogenous and paraphyletic group (to the extent that overall 
similarity mirrors phylogeny). Hellack and Schnell’s analysis also sug¬ 
gested that the genus Saltator, the limits of which have remained essen¬ 
tially unchanged during this century, might require revision. 

With the advent of protein electrophoresis and other molecular-based 
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techniques as tools for determining phylogenetic relationships, we have 
the opportunity to compare the results of genetic data to traditional, 
morphologically based classifications. We chose to focus on the Cardi- 
nalinae because Hellack and Schnell’s classification provided an oppor¬ 
tunity for comparison between a classification based on morphological 
characters with one based on allelic characters. We also sought to clarify 
the placement of the Yellow-shouldered Grosbeak (Caryothraustes hu- 
meralis), because several observers have noted that this species differs in 
some aspects of its natural history from other Caryothraustes species 
(Remsen and Ridgely 1980, Schulenberg et al. 1984, Remsen and Traylor 
1989, Ridgely and Tudor 1989), providing further support for Hellmayr’s 
(1938) concern that this species was not closely related to other Cary¬ 
othraustes. Because skeletal material is now available for several taxa not 
available to Hellack and Schnell, we also include a morphometric analysis 
to examine the phenetic placement of these taxa. 

MATERIALS AND METHODS 

Electrophoresis.— Tissue samples from 27 individuals representing 16 taxa were collected 
between 1981-1985 by personnel of the Museum of Natural Science, Louisiana State Univ. 
Atlapetes brunneinucha (Fringillidae), Catamblyrhynchus diadema (Catamblyrhynchidae), 
and Querula purpurata (Cotingidae) were designated outgroups. The outgroup taxa vary in 
their phylogenetic relationships with the Cardinalinae, and these three outgroups were se¬ 
lected to maximize flexibility in the analyses and to allow comparisons between results 
generated by using different outgroups. Cardinaline taxa were selected to maximize potential 
comparisons with Hellack and Schnell’s study by including at least one taxon from each of 
Hellack and Schnell’s “best phenetic classification” clusters (Table 1). Generally, two in¬ 
dividuals per taxon were examined. Thus, the number of species-level taxa was emphasized, 
rather than number of individuals per taxon, to estimate patterns of genetic variation among 
a variety of taxa. 

Detailed procedures for collection, transport, and storage of tissues followed Johnson et 
al. (1984). Samples consisted of pooled portions of heart, liver, kidney, and pectoral muscle 
(< 1.5 g total) minced and homogenized with an approximately equal volume of deionized 
water, and centrifuged at 12,000 rpm for 30 min. The supernatant was removed and stored 
at — 70°C, and the pellet discarded. Standard horizontal starch-gel electrophoresis (Harris 
and Hopkinson 1976) was performed on the supernatant, using an 11.5% starch-gel and a 
variety of buffer systems (Tris-citrate II, pH 8.0; Lithium hydroxide, pH 8.0; phosphate 
citrate, pH 6.8; and Poulik, pH 8.7; see Selander et al. 1971). 

Loci analyzed were: ADA (E.C. 3.5.4.4), AGPD (1.1.1.8), AK (2.7.4.3), CK (2.7.3.2), 
GDH (1.1.1.47), GOT1 (2.6.1.1), GOT2 (2.6.1.1), ICD1 (1.1.1.42), ICD2 (1.1.1.42), LDH 
(1.1.1.27), PEP-LGG (3.4.11), MDH1 (1.1.1.37), MDH2 (1.1.1.37), MPI (5.3.1.8), NP 
(2.4.2.1), PGI (5.3.1.9), PGM (2.7.5.1), 6PGD (1.1.1.44). Staining techniques and enzyme 
nomenclature follow Harris and Hopkinson (1976). Genotypic data for each individual at 
each locus were entered into Swofford and Selander’s (1981) BIOSYS-1 computer program 
which was used to compute genetic distance coefficients of Nei (1978) and Rogers (1972), 
to derive distance Wagner trees, and to compile associated phenograms based on Rogers’ 
(1972) genetic distance. Strict consensus trees were generated using PAUP (Swofford 1989), 
coding loci as characters and alleles as character states; the stability of the branching diagrams 
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Table 1 

Taxa Studied, Identification Number for PCA Plots, and Sample Sizes for the 

Skeletal and Electrophoretic Analyses 

Num 


Skeletal data 

Electro¬ 

phoresis 

ber 

Species 

N (6) 

N ( 9 ) 

N 

1 

Dickcissel (Spiza americana ) 

1 

1 

2 

2 

Yellow Grosbeak (Pheucticus aureoventris) 

2 

0 

1 

3 

Yellow-bellied Grosbeak (P. chrysogaster) 

0 

1 

0 

4 

Yellow Grosbeak (P. chrysopeplus) 

2 

0 

0 

5 

Rose-breasted Grosbeak (P. ludovicianus) 

6 

3 

2 

6 

Black-headed Grosbeak (P. melanocephalus ) 

2 

2 

0 

7 

Northern Cardinal (Cardinalis cardinalis) 

10 

10 

0 

8 

Pyrrhuloxia (C. sinuatus) 

2 

0 

2 

9 

Yellow-green Grosbeak (Caryothraustes 
canadensis) 

2 

2 

0 

10 

Yellow-shouldered Grosbeak (C. humeral is) 

1 

0 

1 

11 

Crimson-collared Grosbeak (Rhodothraupis 
celaeno) 

0 

2 

0 

12 

Slate-colored Grosbeak (Pitylus grossus) 

4 

1 

1 

13 

Streaked Saltator (Saltator albicollis) 

2 

2 

2 

14 

Black-throated Saltator (S. atricollis) 

2 

0 

0 

15 

Golden-billed Saltator (S. aurantiirostris) 

4 

4 

2 

16 

Masked Saltator (S. cinctus) 

0 

1 

0 

17 

Grayish Saltator (S. coerulescens) 

6 

3 

2 

18 

Buff-throated Saltator (S. maximus ) 

7 

6 

2 

19 

Black-cowled Saltator (5. nigriceps) 

0 

1 

0 

20 

Rufous-bellied Saltator (S. rufiventris) 

2 

0 

0 

21 

Green-winged Saltator (S. similis) 

1 

0 

0 

22 

Lazuli Bunting ( Passerina amoena) 

0 

1 

0 

23 

Ultramarine Grosbeak ( P. brissonii) 

1 

1 

0 

24 

Blue Grosbeak (P. caerulea) 

2 

2 

2 

25 

Painted Bunting (P. ciris) 

2 

2 

0 

26 

Indigo Bunting (P. cyanea ) 

2 

2 

2 

27 

Blue-black Grosbeak (P. cyanoides) 

2 

2 

2 

28 

Indigo Grosbeak (P. glaucocaerulea) 

1 

1 

0 

29 

Orange-breasted Bunting (P. leclancherii) 

1 

0 

0 

30 

Blue Bunting (P. parellina) 

0 

1 

0 

31 

Rose-bellied Bunting (P. rositae) 

1 

0 

0 

32 

Varied Bunting (P. versicolor) 

1 

0 

0 

33 

Rufous-sided Towhee ( Pipilo erythrophthalmus) 

1 

1 

0 

34 

Chipping Sparrow ( Spizella passerina) 

1 

1 

0 

35 

Purple-throated Fruitcrow ( Querula purpurata) 

0 

0 

1 

36 

Chestnut-capped Brush-finch ( Atlapetes 
brunneinucha) 

0 

0 

1 

37 

Plush-capped Finch ( Catamblyrhynchus diadenia) 

0 

0 

1 
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was tested with the bootstrap procedure (Felsenstein 1985). Allelic frequency data generated 
from BIOSYS-1 were entered into the SAS (SAS Institute, Inc. 1985) program “PRIN- 
COMP,” which performed a principal components analysis (PCA) on arcsine square-root 
transformations of allelic frequencies. The transformations are intended to remove the 
dependence of the variance on the mean allelic frequency. Barrowclough and Johnson (1988) 
demonstrated the usefulness of executing a PCA on gene frequency data, showing that, unlike 
morphometric data, the genetic loci analyzed are the products of separate genomic sequences 
and thus are independent of each other. In light of this, PCA will not reduce the dimen¬ 
sionality of the data unless nongenetic factors, such as isolation, gene flow, selection, or 
“phylogenetic relatedness” influence many (or all) of the loci examined. 

Morphometries.—Table 1 lists all species analyzed, the code given to each, the number 
of skeletons, and the composition of samples with regard to sex. The 15 measurements 
selected were those found to be relatively repeatable (in the Fox Sparrow [Passerella iliaca ]) 
by Zink (1983) and were chosen to represent a variety of regions, peripheral and core, of 
the bird skeleton. These measurements, made to the nearest 0.01 mm, are (1) SW = skull 
width, (2) SL = skull length, (3) CO = coracoid length, (4) SC = width of proximal end of 
scapula, (5) ST = sternum length, (6) PS = posterior synsacrum length, (7) SY = greatest 
width of synsacrum, (8) FP = width of proximal end of femur, (9) FD = width of distal end 
of femur, (10) FL = femur length, (11) TB = tibiotarsus length, (12) HT = head of trochanter 
(humerus), (13) HL = humerus length, (14) UL = ulna length, and (15) UP = width of 
proximal end of ulna. The majority of these measurements are described and illustrated by 
Robins and Schnell (1971). 

These linear measurements were analyzed using univariate and multivariate techniques. 
To determine the feasibility of pooling sexes of the same species in the analyses, we quantified 
sexual dimorphism in the two species for which we had the largest samples, the Northern 
Cardinal ( Cardinalis cardinalis) and ( Saltator aurantiirostris). Analysis of variance (ANO- 
VA) using SAS (SAS Institute, Inc. 1985) was used to analyze means of male, female, and 
pooled samples. 

Taxonomic distances, calculated from standardized data, were used to construct UPGMA 
phenograms using the computer program “NTSYS” (Rohlf et al. 1974). Again, the analyses 
were performed on the three different data sets (male, female, and pooled). For each tree, 
a cophenetic correlation coefficient was calculated to indicate the goodness of fit of each 
tree to the original distance matrix (Rohlf et al. 1974). 

Four taxa were represented by specimens lacking the rhampotheca: Passerina rositae, P. 
glaucocaerulea, P. versicolor, and Rhodothraupis celaeno. All analyses were performed with 
and without these taxa. In addition, the only character that would have been affected by 
the absence of the rhampotheca, skull length, was dropped in separate analyses. Results of 
these analyses failed to reveal any significant effect on the placement of the four taxa in 
question. However, omission of skull length did affect the placement of the other taxa. For 
this reason skull length was retained in the remaining analyses. 

Two species of Emberizinae, a different subfamily within the Emberizidae, were included 
as “outgroups”: Chipping Sparrow ( Spizella passerina) and Rufous-sided Towhee ( Pipilo 
erythrophthalmus). The Chipping Sparrow is similar in body size and bill shape to many 
cardinaline buntings, and the towhee is, likewise, superficially similar to many cardinaline 
grosbeaks and saltators. 

Because taxa of the Cardinalinae range from 20 to 80 g in body mass, we attempted to 
correct for the effect of body size on the raw data. Body mass itself was not used because 
of its inherent variability in small samples and because it was not available for several taxa. 
Wiedenfeld (1978) found that for the Tyrannidae the length of the humerus was strongly 
correlated (r = 0.98) with cube root of body weight, and that this was the highest correlation 


Tamplin et al. • BIOCHEMICAL RELATIONSHIPS OF CARDINALINAE 97 


Table 2 

Distribution of Electromorphs (Denoted by Lower Case Letters) in 
Cardinalinae and One Outgroup Taxon ( Qvervla purpurata) 

13 











Locus J 









Species 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

Spiza a me ri can a 
Pheuticus aurco- 

ee 

e 

a 

a 

a 

ac 

a 

d 

b 

a 

e 

a 

a 

a 

df 

a 

b 

d 

ventris 

e 

a 

b 

a 

a 

e 

a 

a 

b 

a 

e 

a 

a 

a 

d 

a 

b 

e 

P. ludovicianus 
Cardinal is sin- 

c 

a 

b 

a 

a 

e 

a 

b 

b 

a 

e 

a 

a 

a 

d 

a 

b 

be 

uatus 

Caryothraustcs 

e 

a 

b 

a 

a 

e 

a 

b 

b 

a 

e 

a 

a 

a 

e 

a 

b 

c 

humeralis 

d 

e 

a 

a 

a 

b 

a 

b 

a 

a 

d 

a 

a 

a 

a 

a 

e 

e 

Pitylus grossus 

e 

b 

b 

a 

a 

e 

a 

d 

b 

a 

e 

a 

a 

a 

bd 

a 

ae 

b 

Saltator albicollis 

e 

b 

b 

a 

a 

c 

a 

ed 

b 

a 

ae 

a 

a 

a 

e 

a 

e 

b 

S. anrantiirostris 

e 

b 

b 

a 

a 

e 

a 

b 

b 

a 

e 

a 

a 

a 

e 

a 

e 

b 

S. coendescens 

e 

b 

b 

a 

a 

e 

a 

d 

b 

a 

e 

a 

a 

a 

e 

a 

c 

ab 

S. maxinms 
Passerina 

e 

b 

b 

a 

a 

e 

a 

e 

b 

a 

ce 

a 

a 

a 

ce 

a 

be 

b 

caerulea 

e 

b 

b 

a 

a 

c 

a 

b 

b 

a 

e 

a 

a 

a 

d 

a 

b 

e 

P. cyanea 

b 

b 

b 

a 

a 

a 

a 

b 

b 

a 

c 

a 

a 

a 

d 

a 

b 

c 

P. cyanoides 
Querula pur- 

e 

b 

b 

a 

a 

e 

a 

b 

a 

a 

ae 

a 

a 

ab 

d 

a 

b 

be 

purata 

a 

d 

c 

a 

a 

d 

b 

e 

c 

a 

b 

a 

a 

e 

d 

a 

b 

e 


3 Loci in order are: ADA, AGDP. AK. CK, GDH, GOT I, GOT2, ICDI, ICD2, LDH, PEP-LGG, MDH1, MDH2, MPI. 
NP, PGI, PGM, 6PGD. 


with weight among 58 skeletal eharaeters. By using the cube root, weight was transformed 
to a linear variable comparable to humerus length (Wiedenfeld 1978). To determine if the 
relationship between humerus length and cube root of body mass holds for the Cardinalinae, 
we calculated the regression coefficient of humerus length on the cube root of weight. This 
regression indicated a strong relationship (r = 0.91). Using humerus length to estimate body 
size, linear regressions were performed with humerus length as the independent variable 
and all other characters as dependent variables. Following Reist's (1986) method for size 
correction, residuals were used as raw data in the same analyses described above. 

RESULTS 

Electrophoresis. —Eighteen loci were identified for each of the individ¬ 
uals sampled, of which five were monomorphic, six showed few (2-3) 
alleles, and seven were highly variable. The distribution of alleles for 18 
loci among the Cardinalinae taxa and one outgroup taxon is presented in 
Table 2. Nei’s (1978) genetic distance values within the Cardinalinae range 
from 0.019 to 0.732 (Table 3). Rogers’ (1972) genetic distance values 
within the Cardinalinae range from 0.056 to 0.529 (Table 3). The average 







Table 3 

Matrix of Genetic Distance Coefficients 1 
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- Spiza artiericana 

- Passerina cyanea 

- - Passerina (Cyanocompsa) cyanoides 

- Pheucticus ludovicianus 

^ ■ ■ - Passerina (Guiraca) caerulea 

r “ Pheucticus aureoventris 

- Cardinalis sinuatus 

- Pitylus grossus 

" Sal tat or coerulescens 
- Saltator albicollis 

__ $ a \ tator aurantiirostris 

- Saltator maxim us 

- Caryothraustes humeralis 

-—- Querula purpurata* 

\ _I_1_I_1-1 

0.60 0.48 0.36 0.24 0.12 0.00 

Genetic distance 

Fig. 1. UPGMA phenogram derived from Rogers’ (1972) genetic distance matrix. The 
cophenetic correlation coefficient is 0.964. Asterisk (*) denotes outgroup. 


intrageneric genetic distance values (D) are: Saltator = 0.066 (N = 4); 
Passerina = 0.013 (N = 3); and Pheucticus = 0.119 (N = 2). The average 
genetic distance between members of the subfamily is 0.324; however, 
computing D for the subfamily excluding Caryothraustes humeralis re¬ 
duces the value to 0.269. 

A UPGMA phenogram (Fig. 1) reveals relationships that are funda¬ 
mentally similar to Hellack and Schnell’s (1977) “best phenetic classifi¬ 
cation.” Querula purpurata was chosen as the outgroup in the distance 
analysis because it displayed the most unique alleles among the potential 
outgroups and it clarified the branching diagrams. Most noteworthy is 
the removal of Caryothraustes humeralis from the Cardinalinae lineage. 
The other major difference is the placement of Cardinalis sinuatus, which 
although phenetically linked to the Saltator complex, clusters genetically 
with the Passerina-Pheucticus complex. 

A 50% majority rule consensus tree, 49 steps long, was produced from 
172 most parsimonious trees (Fig. 2). Consensus trees were generated 
using Atlapetes brunneinucha, Catamblyrhynchus diadema, and Querula 
purpurata as outgroups. These trees show similarities to each other as 
well as with the UPGMA phenogram (Fig. 1). The major differences 
between the consensus tree with Q. purpurata as the outgroup and Fig. 1 
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consistency index is 0.677, excluding uninformative characters. The numbers indicated on 
the branches correspond to the percentage of 172 most parsimonious trees which support 
the placement of that branch. Asterisk (*) denotes outgroup. 


were the placement of Passerina ( Cyanocompsa) cyanoides, which clusters 
with Caryothraustes humeralis due to the presence of a synapomorphy at 
the ICD-2 locus. Caryothraustes humeralis, however, shows eight auta- 
pomorphies, which accounts for its removal from the subfamily in the 
distance analysis. In addition, the placement of Spiza americana within 
the Pitylus-Saltator complex contrasts with its outlying position in Fig. 
1. A strict consensus tree, 41 steps long, with Atlapetes brunneinucha as 
the outgroup is nearly identical to Fig. 2. The consensus tree with C. 
diadema as the outgroup is identical to Fig. 2, with the exception of Spiza 
americana, which is the outlier to all Cardinalines in both this tree as 
well as in Fig. 1 (with the exception of Caryothraustes humeralis). A strict 
consensus tree 50 steps long was generated using both A. brunneinucha 
and Q. purpurata as outgroups. This tree is identical (with the addition 
of Q. purpurata) to the consensus tree with Atlapetes as the only outgroup. 
Using different outgroups and combinations of outgroups had little effect 
on the results of the cladistic analyses. Cladistic analyses of allozyme data 
are compromised, however, because they do not consider gene frequency 
information (Swofford and Berlocher 1987). Bootstrapping (Felsenstein 
1985) the consensus trees collapsed all nodes of the branching diagram, 
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Fig. 3. Distribution of 14 taxa across Principal Components 1 , 2, and 3 based on 
electrophoretic data. Species numbers are: (1) Spiza americana\ (2) Pheucticus aureoventris\ 
(5) P. ludovicianus ; (8) Cardinalis sinuatus ; (10) Caryothraustes humeralis ; (12) Pitylus 
grossus\ (13) Saltator albicollis\ (15) S. aurantiirostris; (17) 5". coerulescens\ (18) S', maximus: 
(24) Passerina caerulea\ (26) P. cyanea\ (27) P. cyanoides; (35) Querula purpurata. 


except the branches associated with the outgroup(s), resulting in one large 
Cardinalinae cluster. 

Principal component 1, which explains 28% of the total variation, sep¬ 
arates the Saltator-Pitylus complex from the Pheucticus-Passerina cluster 
and demonstrates the distinctiveness of the outgroup Querula purpurata 
(Fig. 3). Principal component 2 explains 21% of the total variance and 
serves chiefly to isolate the outgroup from the Cardinalinae, although the 
Passerina-Pheucticus grouping appears well-removed from the remaining 
taxa. Principal component 3 (16% of the total variance) shows marked 
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disjunction between Caryothraustes and other Cardinalinae, and further 
contributes to isolating the outgroup. Thus, in the plot of the first three 
components, which accounts for nearly 65% of the total variance, Car¬ 
yothraustes humeralis and Querula purpurata are clearly distinct from the 
other Cardinalinae, as well as from each other. Loci showing substantial 
contributions of the original variables to the PC axes are: PGM, ADA, 
6GPD, and NP (PC 1); ICD1, GOT1, GOT2, MPI, AK, and AGPD (PC 
2); and ICD2, PGM, NP, PEP-LGG and AK (PC 3). 

Morphometries. — ANOVA revealed pronounced sexual dimorphism in 
skeletal characters. In Cardinalis cardinalis, 10 of 15 characters differed 
significantly (P < 0.05), with the males being larger. In Saltator auran- 
tiirostris, 11 characters differed significantly, with the females having larg¬ 
er values. Differences of this magnitude, producing bimodality of the 
sample means, indicate that the sexes should not be pooled (Wallace 
1984). Because of limited sample sizes, Hellack and Schnell (1977) were 
forced to pool sexes in several species. 

The amount of variation described by PC 1 in the combined male and 
female PCA uncorrected for size ranged from 88.6% in females to 97.4% 
in males. PC 2 described from 2.0% to 5.4% and PC 3, from 0.8% to 
2.7%. Character loadings were equal and positive on PC 1; this suggests 
that PC 1 is strongly related to body size (Wiley 1981, Freeman and 
Jackson 1990). Examination of PCA plots (Figs. 4 and 5) further indicates 
that the effect of size on this analysis is great. Smaller species in the 
expanded genus Passerina are found on the low end of PC 1, whereas the 
largest saltators and grosbeaks are on the higher end. 

UPGMA phenograms without size correction (Figs. 6 and 7) also show 
the effects of body size. The smaller Spizella passerina tends to cluster 
with the smaller cardinalines and Pipilo erythrophthalmus clusters with 
some saltators and other large species. 

In analyses with the data standardized for body size (Fig. 8), PC 1 
describes from 28% to 33% of the total variation, PC 2 ranges from 20% 
to 33%, and PC 3 describes 13% in all three plots. The total variation 
described by the plots is 61% for the combined samples, 63% for the male 
samples, and 69% for the female samples. 

After correction for body size, the two outgroups are well separated 
from the main group, although Saltator nigriceps is found near Pipilo 
erythrophthalmus in both the plots of females and the combined samples. 
Spiza americana is also found outside the subfamily. Pheucticus melan- 
ocephalus and P. ludovicianus are close together as would be expected of 
two allospecies that hybridize (West 1962). However, the allospecies P. 
chrysogaster and P. chrysopeplus did not cluster together, presumably 
because these taxa are only represented by female and male samples 
respectively. The Cardinalis samples, however, do not cluster together in 


Tamplin el al. • BIOCHEMICAL RELATIONSHIPS OF CARDINALINAE 103 


Pc 3 



Fig. 4. Distribution of 28 taxa across Principal Components 1, 2, and 3 based on 
uncorrected morphometric data of male specimens only. Species numbers are: (1) Spiza 
americana: (2) Pheucticus aureoventris: (4) P. chrysopeplus: (5) P. ludovicianus: (6) P. me- 
lanocephalus: (7) Cardinalis cardinally (8) C. sinuatus: (9) Caryothrausles canadensis , (10) 
C. humeralis: (12) Pitylus grossus; (13) Saltator albicollis: (14) S. atricollis ; (15) S. auran- 
tiirostris ; (17) S. coeridescens ; (18) S. maximus: (20) S. rufventris: (21) S’, si mi lis: (23) 
Passerina brissonii: (24) P. caerulea ; (25) /\ c/m; (26) / J . cyanea ; (27) P. cyanoides: (28) P. 
glaucocaerulea ; (29) P. leclancherii ; (31) P. rositae ; (32) P. versicolor ; (33) P/p/Vo erythro- 
phthalnius ; and (34) Spizella passerina. Bold circles indicate taxa used in the electrophoretic 
analysis. 

the plot of male-only samples. Caryothrausles humeralis is found within 
the main group but not near its supposed congener C. canadensis. 

The combined male and female UPGMA phenogram (Fig. 9) corrected 
for size gives results similar to those of the size corrected PCA. Both 
outgroups cluster outside most Cardinalinae, although S. nigriceps again 
is also found outside the subfamily. Caryothrausles humeralis is in the 
main group, but again apart from C. canadensis. The cardinals remain a 
distinct cluster. Pheucticus melanocephalus and P. ludovicianus again clus- 
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Fig. 5. Distribution of 24 taxa across Principal Components 1, 2, and 3 based on 
uncorrected morphometric data of female specimens only. Species numbers are: (1) Spiza 
americana\ (3) Pheucticus chrysogaster, (5) P. ludovicianus\ (6) P. melanocephalus\ (7) Car- 
dinalis cardinalis ; (9) Caryothraustes canadensis ; (11) Rhodothraupis celaeno ; (12) Pity/us 
grossus ; (13) Salta tor albicollis; (15) S. aurantiirostris\ (16) S. cinctus\ (17) S'. coerulescens\ 
(18) S', maximum ; (19) S’. nigriceps\ (22) Passerina amoena; (23) P. brissonr, (24) P. caerulea-, 
(25) P. ciris; (26) P. cyanea\ (27) P. cyanoides\ (28) P. glaucocaerulea; (30) P. parellina\ (33) 
Pipilo erythrophthalmus\ and (34) Spizella passerina. Bold circles indicate taxa used in the 
electrophoretic analysis. 


ter closely. Species from two large genera, Passerina and Saltator, show 
little tendency to cluster together. 

Results from the analyses without size-correction resemble both the 
Hellack-Schnell and the biochemical classifications. However, our results 
indicate that both Passerina ( Guiraca ) caerulea and P. ( Cyanocompsa) 
cyanoides lie outside the main group of cardinalines along with Cary¬ 
othraustes, whereas Hellack and Schnell placed P. caerulea outside Pas- 






























Tamplin et al. • BIOCHEMICAL RELATIONSHIPS OF CARDINALINAE 105 


rC 


Spiza americana 
Cardinalis sinuatus 
Caryothraustes humeralis 
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6.013 5.213 4.413 3.613 2.813 2.013 1.213 0.413 

Taxonomic distance 

Fig. 6. UPGMA phenogram derived from uncorrected morphometric data based on 
male specimens only. The cophenctic correlation coefficient is 0.832. Asterisk (*) denotes 
outgroups. 


serina. In addition, P. brissonii clusters with P. caerulea and P. cyanoides 
in the separate phenograms for the sexes. Hellack and Schnell’s results 
also suggest that Cardinalis and Pheucticus contain species that morpho¬ 
logically are more similar to each other than to species in other genera, 
whereas our results do not indicate this. The fragmented nature of the 
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Taxonomic distance 

Fig. 7. UPGMA phenogram derived from uncorrected morphometric data based on 
female specimens only. The cophenetic correlation coefficient is 0.759. Asterisk (*) denotes 
outgroups. 


genus Saltator results from the inclusion of many taxa not included in 
Hellack and Schnell’s analysis of skeletal characters. The clustering of 
Pitylus grossus and Rhodothraupis celaeno within the saltators agrees with 
Hellack and Schnell’s (1977) results. The results of the size-corrected 
phenograms agree with our biochemical results on the distinctiveness of 
Spiza americana. However, the separation of two pairs of hybridizing 
allospecies (Passerina cyanea and P. amoena; Pheucticus ludovicianus and 
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Fig. 8. Distribution of 34 taxa across Principal Components 1, 2, and 3 based on size- 
corrected morphological data of both males and females. Species numbers are: (1) Spiza 
americana: (2) Pheucticus aureoventris: (3) Phcucticus chrysogaster ; (4) P. chrysopeplus: (5) 
P. ludovicianus ; (6) P. melanocephalus: (7) Cardinalis cardinalis: (8) C. sinuatus: (9) Car- 
yothraustes canadensis’, (10) C. humeralis ; (11) Rhodothraupis celaeno: (12) Pitylus grossus; 
(13) Saltator albicollis ; (14) S. atricollis\ (15) S. aurantiirostris: (16) S. cinctus; (17) S. 
cocndcscens ; (18) 5. maximus\ (19) S. nigriccps ; (20) S. rufiventris\ (21) S. si mi/is: (22) 
Passerina amoena ; (23) P. brissoni ; (24) P. caeru/ca ; (25) P. c/m; (26) P. cvanea: (27) P. 
cyanoidcs: (28) P. glaucocaerulea ; (29) P. lec/anchcrii ; (30) P. parellina: (31) P. rositae ; (32) 
P. versicolor ; (33) Pipilo erythrophthalmus ; and (34) Spizella passerina. Bold circles indicate 
taxa used in the electrophoretic analysis. 

P. melanocephalus) indicates problems with phylogenetic interpretations 
of these results. 


DISCUSSION 

The results of our morphometric analysis present a problem to system- 
atists who would use morphometries to infer phylogeny. If correction for 
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Fig. 9. UPGMA phenogram derived from size-corrected morphological data of both 
sexes. The cophenetic correlation coefficient is 0.789. Asterisk (*) denotes outgroups. 
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size is not made, then it is impossible to distinguish outgroups. However, 
if correction for size is made, then the resulting phenogram makes little 
sense. For example, the phenogram places two hybridizing allospecies, 
Passerina amoena and P. cyanea (Sibley and Short 1959), in different 
clusters. Phenetic analyses of morphometric data may not be able to 
separate outgroups when variability in size is large among the taxa under 
investigation and when shape differences are small, as in the Cardinalinae, 
which seem to have undergone little divergence in shape. On the other 
hand, our results suggest that size itself has a phylogenetic basis within 
the Cardinalinae. Although we cannot rule out the possibility that mea¬ 
surement of additional characters would have provided more resolution, 
most morphometric analyses of birds have shown that most skeletal char¬ 
acters, particularly in the trunk skeleton, are highly autocorrelated and 
contribute little additional information. 

Relationships of Passerina, Guiraca, and Cyanocompsa.— Hellmayr 
(1938) and Hellack and Schnell (1977) favored retention of all three gen¬ 
era, whereas Paynter (1970) merged Guiraca and Cyanocompsa in Pas¬ 
serina without comment. Because biochemical data suggest that Passerina 
{Guiraca) caerulea might be more closely related to Pheucticus ludovici- 
anus than to either Passerina cyanea or Passerina ( Cyanocompsa) cy- 
anoides, our results support the classifications of Hellmayr (1938), Hellack 
and Schnell (1977), the AOU (1983), Ridgely and Tudor (1989), and 
Sibley and Monroe (1990), in retaining the genera Guiraca and Cyano¬ 
compsa. Unfortunately, tissue samples for the Blue Bunting {Passerina 
[ Cyanocompsa ] parellina), the Indigo Grosbeak {Passerina [Cyanoloxia] 
glaucocaerulea), and the Blue Finch {Passerina [Porphyrospiza] caerules- 
cens ) were not available. Ridgely and Tudor (1989) removed the latter 
species from the Cardinalinae, a change supported by Bates et al. (in 
press). 

Monophyly of Pheucticus.— Hellack and Schnell (1977) found that the 
genus Pheucticus as viewed by Paynter (1970) consisted of two rather 
different groups, one of which was more closely related to Guiraca (see 
above). Hellmayr (1938) had previously separated these two groups into 
two genera, retaining Hedymeles for ludovicianus and melanocephalus. 
Our results, both biochemical and morphometric, support the classifi¬ 
cations of Hellmayr and Hellack and Schnell. Paynter’s broader Pheucticus 
may be paraphyletic, with Guiraca seemingly more closely related to 
“ Pheucticus ” ludovicianus and Cardinalis cardinalis more closely related 
to P. aureoventris. However, the difference in allozymes is so small among 
these taxa that we hesitate to recommend resurrection of Hedymeles for 
ludovicianus and melanocephalus without corroborating data from other 
kinds of analyses. The strong plumage and vocal similarities between 
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ludovicianus and melanocephalus and the aureoventris group (Ridgely and 
Tudor 1989; Remsen, pers. obs.) seem unlikely to be merely shared prim¬ 
itive characters. Likewise, the strong plumage and vocal similarities of 
Guiraca caerulea to the Passerina buntings seem unlikely to be merely 
shared primitive characters. These similarities, obviously, have been re¬ 
sponsible for current taxonomic arrangement. 

Monophyly of Saltator. — Hellack (1976) and Hellack and Schnell (1977) 
concluded that Saltator is not a monophyletic group. Our results also 
support this conclusion but for different reasons. Our biochemical data 
indicate that Pitylus grossus, traditionally treated as a close relative of 
Saltator, is more closely related to S. coerulescens than the latter is to 
three other species of Saltator (Fig. 1). Hellack and Schnell (1977) also 
found that P. grossus was more closely related to some saltators, including 
S. coerulescens, than those saltators were to other Saltator species; their 
branching pattern, however, differed in several ways from that shown in 
Fig. 1. For example, they found that Saltator maximus clustered with S. 
coerulescens, whereas we found that S. maximus was the most divergent 
Saltator analyzed. Because we analyzed only four species within this large 
genus, specific taxonomic recommendations are not warranted other than 
the merger of Pitylus into Saltator. 

Relationships of Spiza americana. —Although current classifications 
(Morony et al. 1975, AOU 1983) follow Paynter (1970) and Hellack and 
Schnell (1977) in placing the Dickcissel in the Cardinalinae, others have 
proposed that its relationships lie elsewhere. Our results indicate that 
Spiza americana is the outgroup to all Cardinalinae analyzed except Car- 
yothraustes humeralis, which is probably not a true cardinaline (see be¬ 
low). This was also the case in our size-corrected morphological data. 
This contrasts with the results of Hellack and Schnell (1977), who found 
that Spiza was the sister taxon to Passerina sensu strictu. Clearly, a more 
extensive analysis is needed to determine the relationships of Spiza. 

Relationships of Caryothraustes humeralis.—The Yellow-shouldered 
Grosbeak (C. humeralis) has been included with the other species in the 
genus Caryothraustes since the classification of Hellmayr (1938). Hellmayr 
noted, however, that C. humeralis “probably deserves generic separation” 
and gave plumage and bill shape features that distinguished it from other 
Caryothraustes. Hellmayr proposed that it was “intermediate” between 
Saltator, the genus in which it was placed by Chapman (1926), and the 
other Caryothraustes. Hellack and Schnell (1977), however, retained hu¬ 
meralis in Caryothraustes on the basis of plumage characters, and their 
analysis suggested that Caryothraustes was most closely related to Cyano- 
compsa. 

Virtually nothing was known of the natural history of C. humeralis until 
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recent fieldwork in western Amazonia revealed that the species was more 
widely distributed than formerly believed and that it was a scarce member 
of mixed-species canopy flocks composed primarily of tanagers (Remsen 
and Ridgely 1980, Schulenberg and Remsen 1982, Cardiff 1983, Schu- 
lenberg et al. 1984). We propose that humeralis might not be related to 
other Caryothraustes species because, in addition to the characters pointed 
out by Hellmayr (1938), humeralis differs from other Caryothraustes 
members in its feeding social system (quiet, inconspicuous, and solitary 
or in pairs in mixed-species flocks versus noisy, single-species flocks). 

Our results suggest that Caryothraustes humeralis is not a member of 
the Cardinalinae. It is the outgroup to all taxa analyzed, including Spiza, 
other than the cotinga Querula purpurata. In terms of Nei’s (1978) genetic 
distance, it is almost as distant from the Cardinalinae as is Querula, a 
member of a different suborder within the Passeriformes. 

In the absence of tissue samples from other Caryothraustes species, the 
relationships of humeralis to its purported congeners cannot be ascer¬ 
tained at present. It is not surprising that C. humeralis clustered with its 
congener, C. canadensis, in the morphological phenograms uncorrected 
for size because overall similarities in size and shape probably contributed 
to the original allocation of humeralis to the genus Caryothraustes. We 
predict that humeralis will be found to be unrelated to other species in 
the genus Caryothraustes. 
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